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We present an exhaustive analysis of scalar, electromagnetic and gravitational perturbations in the 
background of Schwarzchild-de Sitter and Reissner-Nordstrom-de Sitter spacetimes. The field prop- 
agation is considered by means of a semi-analytical (WKB) approach and two numerical schemes: 
the characteristic and general initial value integrations. The results are compared near the ex- 
treme cosmological constant regime, where analytical results are presented. A unifying picture is 
established for the dynamics of different spin fields. 
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I. INTRODUCTION 

Wave propagation around nontrivial solutions of Ein- 
stein equations, black holes in particular, is an active field 
of research (see Q, H, |^ and references therein). The 
perspective of gravitational wave detection in the near 
future and the great development of numerical general 
relativity have increased even further the activity on this 
field. Gravitational waves should be especially strong 
when emitted by black holes. The study of the propaga- 
tion of perturbations around them is, hence, essential to 
provide templates for gravitational wave identification. 
On the other hand, recent astrophysical observations in- 
dicate that the universe is undergoing an accelerated ex- 
pansion phase, suggesting the existence of a small posi- 
tive cosmological constant and that de Sitter (dS) geome- 
try provides a good description of very large scales of the 
universe 0. We notice also that string theory has re- 
cently motivated many works on asymptotically anti-de 
Sitter spacetimes (see, for instance, |EIE0I1|)- 

In this work, we perform an exhaustive investigation of 
scalar, electromagnetic, and gravitational perturbations 
in the background of Schwarzchild~de Sitter (SdS) and 
Reissner-Nordstrom-de Sitter (RNdS) spacetimes. Con- 
trasting with the noncharged case, in the RNdS one the 
electromagnetic and gravitational perturbations are nec- 
essarily coupled. We scan the full range of the cosmo- 
logical constant, from the asymptotic flat case (A — 0) 
up to the critical value of A which characterizes, for the 
noncharged case, the Nariai solution Two different 
numerical methods and a higher order WKB analysis 
are used. The results are compared near the extreme 
A regime, where analytical results can be obtained. 

We recall that for any perturbation in the spacetimes 
we consider that, after the initial transient phase, there 
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are two main contributions to the resulting asymptotic 
wave 

[13 113: initially the so-called quasinormal modes, 
which are suppressed at later time by the tails. The 
first can be understood as candidates for normal modes 
which, however, decay (their energy eigenvalues becomes 
complex), as in the ingenious mechanism first described 
by Gamow in the context of nuclear physics [l^ . After 
the initial transient phase, the properties of the result- 
ing waves are more related to the background spacetime 
rather than to the source itself. 

It is well known that for asymptotically flat back- 
grounds the tails decay according to a power law, whereas 
in a space with a positive cosmological constant the decay 
is exponential. Curiously, £ — modes for scalar fields 
in dS spacetimes, contrasting with the asymptotic flat 
cases, exponentially approach a nonvanishing asymptotic 
value 0,[i3|. We detected, by using a noncharacteris- 
tic numerical integration scheme, a dependence of this 
asymptotic value on the initial velocities. In particular, 
it vanishes for static initial conditions. Our results are in 
perfect agreement with the analytical predictions of ■ 

The semianalytical analyses of this work were per- 
formed by using the higher order WKB method proposed 
bv Schutz and Will [l^, and improved by Iyer and Will 
[la. Il7j . It provides a very accurate and systematic way 
to study black hole quasinormal modes. We apply it to 
the study of various perturbation fields in the nonasymp- 
totically flat dS geometry. Quasinormal modes are also 
calculated according to this approximation, and the re- 
sults are compared to the numerical ones whenever ap- 
propriate, providing a quite complete picture of the ques- 
tion of quasinormal perturbations for dS black holes. 

Concerning the charged case, we analyze in detail the 
wave propagation of the massless scalar field and cou- 
pled electromagnetic and gravitational fields in the RNdS 
spacetime. An important difference concerning the dy- 
namics of the electromagnetic and gravitational fields is 
that there are no pure modes, since both are interre- 
lated. We will show that the direct picture of the evolu- 
tion presents us with perfect agreement of quasinormal 
frequencies with those obtained by using the approxima- 
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tion method suggested in [l6] . One important point 
assessed is the dependence of the fields' decay on the elec- 
tric charge of the black hole, including the asymptotically 
flat limit, for which we expected to find traces of a power 
law tail appearing between the quasinormal modes and 
the exponential tail. 

Two very recent works overlap our analysis presented 
here. A similar WKB approach, presented in was 
used very recently by Zhidenko to study SdS black 
holes, giviiigresults in agreement with ours. Yoshida and 
Futamase |2_(| used a continued fraction numerical code 
to calculate quasinormal mode frequencies, with special 
emphasis on high order modes. Our results are also com- 
patible. Finally, we notice that solutions of the wave 
equation in a nontrivial background have also been used 
to infer intrinsic properties of the spacetime ,2 Lj . 

The paper is organized as follows. Sec. II provides 
theoretical considerations and reviews some well-known 
results that were useful to our work; Sec. Ill briefly 
explains the numerical and semianalytical methods em- 
ployed, followed by Sec. IV, which presents, in detail, our 
results on field dynamics for near extreme SdS and RNdS 
geometries. Sec. V deals with the so-called intermediary 
region, where the geometries are not extreme. Data on 
the SdS limit and on exponential tails are also presented. 
Sec. VI deals with the near asymptotically flat region, 
and Sec. VII presents our conclusions. 



II. METRIC, FIELDS, AND EFFECTIVE 
POTENTIALS 



is the SdS geometry in which we are interested. The 
horizons and rc, with r_|_ < Tc, are denoted the event 
and cosmological horizons, respectively. In this case, the 
constants m and a are related to the roots by 
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a =r^+r^ + r+r^ , 



2ma? = r^rc{r+ + Tc) 



(3) 



(4) 



If 9m^A = 1, the zeros r+ and Tc degenerate into a double 
root. This is the extreme SdS black hole. If m^A > 1, 
there are no real positive zeros, and the metric Q does 
not describe a black hole. 

In both the SdS and RNdS cases, we shall study the 
perturbation flelds in the exterior region, defined as 



r+ = {(i,r,0,, 



< r < rd 



(5) 



In this region T_|_, we define a "tortoise coordinate" x{r) 
in the usual way. 



x{r) 



ln(rc - r) + ln(r - r+) 
- -^ln(r - r_) + -^ln(r - r„) , (6) 



with 



dh{r) 



dr 



(7) 



The metric describing a charged, asymptotically de Sit- 
ter spherical black hole, written in spherical coordinates, 
is given by 



-h{r)dt'^ + h{r)-'^dr^ +r^{d9^ +sm^ed(j)^) , (1) 



where the function h(r) is 



, , , 2m ff^ Ar^ 
h{r) = 1 + ^ - 



(2) 



The integration constants m and q are the black hole 
mass and electric charge, respectively. If the cosmological 
constant is positive, we have the Reissner-Nordstr6m~de 
Sitter metric. In this case, A is usually written as A = 
where the constant a is the "cosmological radius." 

The spacetime causal structure depends strongly on 
the zeros of h{r). Depending on the parameters m, g, and 
a, the function h(r) may have three, two, or even no real 
positive zeros. For the RNdS cases we are interested in, 
h(r) has three simple real, positive roots (rc, r^, and r-), 
and a real and negative root r„ — — (r_ + r+ + Tc). The 
horizons r_, r^, and Tc, with r_ < r-|_ < Tc, are denoted 
Cauchy, event, and cosmological horizons, respectively. 

For the SdS case {q = 0), and assuming m > and 
< 9m^A < 1, the function h(r) has two positive zeros 
r+ and Tc and a negative zero r„ = — (r+ + rc). This 



The constants k_, k+, and the surface gravities 

associated with the Cauchy, event, and cosmological hori- 
zons, respectively. For the SdS case, the term associated 
with the Cauchy horizon is absent. 

Consider now a scalar perturbation field $ obeying the 
massless Klein-Gordon equation 



= 



(8) 



The usual separation of variables in terms of a radial field 
and a spherical harmonic Y£.„i(0, ip), 



(9) 



leads to Schrodinger-type equations in the tortoise coor- 
dinate for each value of £, 



+ 



where the effective potential V'^ is given by 



V^'ir) = h{r) 



£{£+!) _^ 2m 2q^ 



(10) 



(11) 



The situation for higher spin perturbations is quite 
different. In the SdS geometry, in contrast to the case 
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of an electrically charged black hole, it is possible to 
have pure electromagnetic and gravitational perturba- 
tions. For both cases, we have Schrodinger-type effective 
equations. For the first, the effective potential is given 
by HI 



h{r) 



U--{pi-p2)W 



(18) 



F^'(r) = h{r) 



(12) 



Pi 



3m + \/ 9rn^ + Scq^ , 



(19) 



with t > 1. The gravitational perturbation theory for 
the exterior Schwarzschild-de Sitter geometry has been 
developed in 0,01. The potentials for the axial and polar 
modes are, respectively, 



P2 



3m — •\/ 9m^ + Scg^ , 



(20) 



y"^(r) = h{r) 



6m 



rh(r) 1 / 2 \ 

(13) W=—^i2cr + 3m) + -[cr + m+-Ar^], (21) 



VP°{r) = 



2h{r) 



(3m, + cr)^ 
+im^{2,cr - Kr^)] , 



[9m^ + 3c^mr^ + c^(l + c)r^ 

(14) 



with 2c = (^-l)(^+2) and £ > 2. For perturbations with 
£ > 0, we can show explicitly that all the effective po- 
tentials V{x) = V{r{x)) are positive definite. For scalar 
perturbations with £ = 0, however, the effective potential 
has one zero point xq and it is negative for x > xq. 

The perturbation theory for the RNdS geometry has 
been developed in 123]. There are neither purely elec- 
tromagnetic nor gravitational modes. Indeed, we have 
four mixed electromagnetic and gravitational fields, two 
of them called polar fields, and Z2 (since they im- 
part no rotation to the black hole) and two named axial 
fields, and Z2 ■ It is possible to express their dynam- 
ics in four decoupled wave equations, two for the axial 
fields and two for the polar fields. Their deduction can 
be found in |2^ and references therein. Here we just 
show the expressions, which will be useful throughout 
our work. 

The axial perturbations Z^2 ^'^^ governed by wave 
equations which have the same form as Eq. Hl()|l . but 
with effective potentials given by 



Ff - h{r) 



(15) 



h{r) 



4g2 3to + y/9m^+8cq^ 

rp'2, ^3 



(16) 

respectively, with 2c = (£ - 1)(^ + 2). 

The polar perturbations Z^2 ^'^^ subjected to rather 
cumbersome potentials, as we can see below: 



h{r) 



U+-{P1-P2)W 



(17) 



2 \ 2cr^h(r) 
U = (2cr -I- 3m,)T4^ + ( w - cr - m Ar^ , 

3 J uj 

(22) 



3m — 



(23) 



In the limit q ^ 0, the RNdS potentials Vf go into 
the SdS polar and axial potentials, . Therefore, the 
minimum £ for these fields is ^ = 2, while the V-^ fields 
admit ^ = 1 as their minimum £ value, since they become 
electromagnetic perturbations in the limit q ^ 0. 



III. NUMERICAL AND SEMIANALYTICAL 
APPROACHES 

A. Characteristic integration 

In [231 a simple but at the same time very efficient 
way of dealing with two-dimensional d'Alembertians has 
been set up. Along the general lines of the pioneering 
work [2^ , the authors introduced light-cone variables u = 
t — X and V = t -\- x^ in terms of which all the wave 
equations introduced have the same form. We call V 
the generic effective potential and ip^ the generic field, 
and the equations can be written, in terms of the null 
coordinates, as 



92 



dudv 



ipi{u,v) = V{r{u,v))ipi{u,v) 



(24) 



In the characteristic initial value problem, initial data 
are specified on the two null surfaces u = uq and v = vq. 
Since the basic aspects of the field decay are independent 
of the initial conditions (this fact is confirmed by our 
simulations), we use 



ipe{u = uo,v) = cxp 



2cr2 



(25) 
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ipeiu, V vq) = cxp 



2a2 



(26) 



Due to the size of our lattices, the latter constant can be 
set to zero for any practical purpose. 

Since we do not have analytic solutions to the time- 
dependent wave equation with the effective potentials in- 
troduced, one approach is to discretize Eq. (|24|l and then 
implement a finite differencing scheme to solve it numer- 
icallv- One possible discretization, used, for example, in 

smiii, is 



MN) = MW) + ME)-MS) 
MW) + ME) 



A'V{S) 



(27) 



where we have used the definitions for the points N = 
{u + A,v + A), W = {u + A,v), E = {u,v + A), and 
S — {u,v). With the use of expression H27|l . the basic 
algorithm will cover the region of interest in the u — v 
plane, using the value of the field at three points in order 
to compute it at a fourth. 

After the integration is completed, the values 
i^i{um.ax,v) and ^iiu,Vmax) are extracted, where Umax 
(vmax) is the maximum value of u (v) on the numerical 
grid. Taking sufhciently large Umax and Vmax, we have 
good approximations for the wave function at the event 
and cosmological horizons. 



We make an exhaustive analysis of the asymptotic be- 
havior of the solutions of Eq. H1U|I with initial conditions 
of the form 



exp 



^'^(O, x) = A exp 



(X -xq) 
2al 



[x - XiY 



2af 



(31) 



(32) 



The results do not depend on the details of the initial 
conditions. They are compatible with the ones obtained 
by the usual characteristic integration, with the only, and 
significant, exception of the £ = scalar mode. As we 
will see, its asymptotic value depends strongly on the 
initial velocities V'£(0,a;), a behavior already advanced in 
the work fill. 



C. WKB analysis 

Considering the Laplace transform of Eq. H10(l . one 
gets the ordinary differential equation 



dx"^ 



[s' + V{x)\^Jt{x)=Q 



(33) 



B. Noncharacteristic integration 

It is not difficult to set up a numeric algorithm to solve 
Eq. (|10l) with Cauchy data specified on a t constant 
surface. We used a fourth order in x and a second order 
in t scheme (see, for instance, |2fij for an application of 
this algorithm to seismic analysis). The second spatial 
derivative at a point {t,x), up to fourth order, is given 
by 

^ [%ljt{t,x + 2Ax) ~l&^t{t,x + Ax) 



Ax) 



-I- 30V'^(i, x) — 16V'£(i, X 
+'ipt{t,x-2Ax)] , (28) 

while the second time derivative up to second order is 

'ipt{t + At, x) - 2ijji{t, x) +i}ji{t- At, x) 



tpi{t,x) 



At2 



(29) 

Given ^^(i = tQ,x) and ipi = {t = to — At,x) 
(or il;i(t — ta,x)), we can use the discretization in 
Eqs. 1211) and it?^ to solve Eq. ((Tn|l and calculate 
tpeit — to + At, x). This is the basic algorithm. At each 
interaction, one can control the error by using the invari- 
ant integral (the wave energy) associated with Eq. IjlUfl 



{■^',{t,x))^ + {i;e{t,x)y + V{x)Mt,x) 



dx . 
(30) 



One finds that there is a discrete set of possible values 
for s such that the function ipi, the Laplace-transformed 
field, satisfies both boundary conditions. 



hm il^e e^^ — 1 , 



X — > — oo 



lim ipi 

x — >+oo 



(34) 



(35) 



By making the formal replacement s = ioj , we have the 
usual quasinormal mode boundary conditions. The fre- 
quencies u) (or s) are called quasinormal frequencies. 

The semianalytic approach used in this work 0, 
is a very efficient algorithm to calculate the quasinormal 
frequencies, which have been applied in a variety of sit- 
uations [13 • With this method, the quasinormal modes 
are given by 



= {Vo + P)-i 



-2V, 



(2) 



1/2 



(1 + g) (36) 



where the quantities P and Q are determined using 



^^8 



Vr 



(4) 



Vr 



(2) 



4+" j-288 



(3) 



Vr 



(2) 



7 + 60a^ 



(37) 
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Q = 



-2Vr 



(2) 



6912 



Vr 



(3) 



Vr 



(2) 



(77+188.^)-^ 



t/(3)2t.(4) 



(2) 



(51 + lOOa^ 



2304 
1 

~ 288 



Vr 



(4) 



(2) 
) 

(6) 



(67 + 68a^ 



(5 + Aa^ 



1 

288 



(3)t^(5) 



vrvr 



Vr 



(2)2 



[l9 + 28a^) 



(38) 



In Eqs. a — n + 1/2 and the superscript 

(«) denotes differentiation, with respect to x, of the po- 
tential V{x). The potential and its derivatives are then 
calculated at the point xq, where V{x) is an extremum. 
The integer n labels the modes 



0,1,2,..., 

-1,-2,-3, 



Re(a;n) > , 
Re(w„) < . 



(39) 



IV. NEAR EXTREME LIMIT 
A. Schwarzschild— de Sitter black hole 

To characterize the near extreme limit of the 
Schwarzschild-de Sitter geometry, it is convenient to de- 
fine the dimensionless parameter S: 



(40) 



The limit < 5 <C 1 is the near extreme limit, where 
the horizons are distinct, but very close. In this regime, 
analytical expressions for the frequencies have been cal- 
culated [23,123. For the scalar and electromagnetic fields, 
the quasinormal frequencies are 



- - 3m' 
3 



1/2 



1/2 




(41) 



For the axial and polar gravitational fields, the frequen- 
cies are given by 



^ o 2 A 2 
3m A 

3 



1/2 



(^ + 2)(^-l)-- 



1/2 





(42) 



1500 



Figure 1: Decay of the scalar and electromagnetic fields with 
I — 1, and of the scalar, axial, and polar gravitational fields, 
with £ — 2, with the SdS geometry approaching the near 
extreme limit. The parameters for the geometry are m — 1.0 
and 5 = 0.01, 0.1, 0.3. 



Direct calculation of the wave functions confirms that, 
in the near extreme limit, their dynamics is simple, with 
the late-time decay of the fields being dominated by 
quasinormal modes. All the types of perturbation tend 
to coincide near the extreme limit. In addition, as we ap- 
proach the extreme limit, the oscillation period increases 
and the exponential decay rate decreases. These conclu- 
sions, illustrated in Fig. [J fo r £ = 1,2, are consistent 
with the ones presented in |27ll2^ . 

By using a nonlinear fitting based on a analysis, it 
is possible to estimate the real and imaginary parts of the 
n = quasinormal mode. These results can be compared 
with the analytical expressions in the near extreme cases. 
In Fig. 121 we analyze the dependence of the frequencies 
on £. The accordance between analytic and numerical 
data is extremely good. 



Reissner— Nordstrom— de Sitter black hole 



They can be compared with the numerical and semiana- 
lytic methods presented in the previous section. 



For the RNdS case, the near extreme limit corresponds 
to the situation where the event and cosmological hori- 
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Figure 2: Curve of Re(a;o)^ with £{£ + 1), for the scalar, 
electromagnetic and gravitational fields, in the near extreme 
limit. The dots are the numerical results and the solid lines 
represent the analytical results. The parameters for the ge- 
ometry are m = 1.0 and S = 0.01. 



zons are very close to each other. It is natural to define 
the dimensionless parameter S as 



(Zsc), two axial fields {Z-^ , Z2), and two polar fields 
Z+). For each one, we have a different potential. 
In the near extreme limit, we have 



Vix) = nir+)h{x) + 0{S) 



Vo 



cosh {k+x) 



-OiS) . (45) 



The constant Vq in the scalar case is denoted by VJ"^, in 
the axial cases by {Vq^ , Vq }^ Siiid in the polar cases by 
{yo^^yo^J^ The fore going expression is a Poschl- Teller 
potential |29l |. 

For the scalar field, Vq'^ has been calculated in psf . 
and is given by 



Vq — 



2r\ 



(46) 



We proceed to the analysis of the coupled electromag- 
netic and gravitational fields. We take the analytical ex- 
pressions for all potentials ^'^'^ to the near 
extreme limit. For the two axial potentials, we have 

^0^' - + l)rl + - r^S,] , (47) 



S = 



rc - r+ 



(43) 



where < 6 ^ 1. In this limit the dynamics can be 
analytically characterized, as has been analyzed in |27j . 
More general settings, including RNdS geometries, were 
explored in p^ . 

The function h{x) = h{r{x)) can be analytically calcu- 
lated j27t^], with the result 



h{x) 



(44) 



2 cosh^(K+x) 

We have five different fields at hand: the scalar field 



where 



^-r^-^[iii+l)rl+Aq'-r+S,] , (48) 



S*! = 3m - ^/9m^ + A{e + 2){e-l)q^ , (49) 



5*2 = 3m -I- y/9m'^ + 4:{£ + 2){£ -l)q^ 



(50) 



We can now turn to the two polar potentials. The 
constants Vq'^''' and Vq~^ are given by 



Vq — 



fi+irc - r+) 



^2cr^ -|- 3m,r-)_ -|- r+ -^/Qm^ + Scg^^ ^cr-|_ + m + —j-^^ 



2r\ 



3m 



c 



(51) 



K+irc ~ r+) 



2cr\_ + 3m,r+ — 



2Arl 



2r\. 



cr+ 



3m 



c 



(52) 



with C = 2mr+ - 2q^ - 2Ar4 /3 and 2c = (£ + 2){l - 1). 
The quasinormal modes associated with the Poschl- Teller 
potential have been extensively studied [s^, The 



frequencies w„ are given by 
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Figure 3: Near extreme fundamental frequencies for differ- 
ent values of the charge q, for the lower multipole mode of 
the scalar field. Analytical values are represented by straight 
lines, and numerical values appear as dots. The parameters 
for the geometry are m = 1.0 and S = 10"'^. 




Figure 4: Near extreme fundamental frequencies for different 
values of the charge q, for the lower multipole mode of the 
polar fields. On the left are the data for and, on the left, 
for Z2 ■ The parameters for the geometry are m = 1.0 and 
5= 10-3. 



with n G {0, 1, . . .} labeling the modes. Using expression 
(|53|l and the expression for Vq, the frequencies can be 
easily calculated. 

We can also use the numerical method to analyze the 
field decay in the near extreme limit. Using a nonlinear 
fitting based on analysis for the wave functions, we 
can estimate the quasinormal frequencies. These results 
can be compared with the analytical expressions in the 
near extreme cases. The accordance between both sets 
of results is extremely good. We illustrate this point in 
Figs. HEl 

Direct calculation of the wave functions confirms that, 
in the near extreme limit, the dynamics of the fields is 
simple, with the late-time decay being completely domi- 
nated by quasinormal modes. 



V. INTERMEDIARY REGION IN PARAMETER 
SPACE 



z'l z". 




Figure 5: Near extreme fundamental frequencies for different 
values of the charge q, for the lower multipole mode of the 
axial fields. The data for appear on the left, and those 
for Z2 , on the right. The parameters for the geometry are 
m = 1.0 and 5 = 10"^ 



A. Schwarzschild— de Sitter black hole 

1. Scalar field with £ = 

Only scalar perturbations can have zero total angular 
momentum. Solutions of Eq. (|10|l with i = lead to 
a constant tail, as already shown in 0, This is 

confirmed in Fig. The novelty here is the dependence 
of the asymptotic value on the ipiiOjX) initial condition. 
Figure El reveals the appearance of the constant value (j)Q 
for large t, and its dependence on ipi{0,x). Note that (po 
falls below 10"'' for ipi{Q,x) = 0. These results are in 
accordance with the analytical predictions of [l^ , which 



give 

^(oo,r) = 4 r^(0,s)s-^. (54) 

fc Jo V^) 



2. Fields with l> 

We can have scalar and vector fields with angular mo- 
mentum £ = 1, and with £ > 1, it is possible to introduce 
gravitational fields also. Their behavior is described in 
general by three phases. The first corresponds to the 
quasinormal modes generated from the presence of the 
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1=2 




Figure 6: Asymptotic solutions <j>{x,t) obtained by nonchar- 
acteristic numerical integration with ^ = 0, A = 10~* and 
m = 1.0. The curves correspond to diflterent values of A in 
the initial condition l|82|l . 



black hole itself. A little later there is a region of power 
law decay, which continues indefinitely in an asymptotic 
fiat space. In the presence of a positive cosmological con- 
stant, however, an exponential decay takes over in the 
latest period. 

As the separation of the horizons increases, the quasi- 
normal frequencies deviate from those predicted by ex- 
pressions and (021) • In Fig- HI this is illustrated 
for ^ = 1,2. Some qualitatively different effects show up 
when we turn away from the near extreme limit. For 
a small cosmological constant the asymptotic behavior 
is dominated by an exponentially decaying mode rather 
than by a quasinormal mode, for all perturbations con- 
sidered. 

It is interesting to compare the values obtained for the 
fundamental modes using the numerical and semianalytic 
methods. We find that the agreement between them is 
good, for the whole range of A. The difference is smaller 
for the first values of £. This is expected, since the nu- 
merical calculations work better in this region. In Table 
U we illustrate these observations with a few values of 
A. It is important to mention that quasinormal frequen- 
cies for the SdS black hole were already calculated in a 
recent paper by applying a variation of the WKB 
method used here There are earlier papers calcu- 

lating quasinormal modes in this geometry, for example, 
112]. 

The first higher n modes cannot be obtained from the 
numerical solution, but can be calculated by the semian- 
alytical method. As the cosmological constant decreases, 
the real and imaginary parts of the frequencies increase, 
up to the limit where the geometry is asymptotically fiat. 
The behavior of the modes is illustrated in Fig. |H1 The 
behavior of the electromagnetic field is similar. 

A analysis of the data presented in Figure El shows 



Figure 7: Graphs of the real and imaginary part of the fun- 
damental frequencies (n — 0) with k+. The dotted lines are 
the near extreme results, the dots are the numerical results 
and the continuous curves are the semianalytic results. In the 
graphs, m = 1.0. 
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Figure 8: Quasinormal modes of the scalar, axial and polar 
gravitational fields, for higher modes. The parameter for the 
curves are m = 1.0 and A — 10"'^. 



that the massless scalar, electromagnetic, and gravita- 
tional perturbations in SdS geometry behave as 



ipl" - e ''^^p* with t oo 



re 



with t oo 



^p* with i — > oo , 



^ e-'^-p* with t 00 



(55) 



(56) 



(57) 



(58) 



for t sufficiently large. At the event and the cosmological 
horizons, t is substituted, respectively, by v and u. 
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Table I: Fundamental frequencies for the scalar field in SdS, obtained using the numerical and semianalytical methods. In this 
table, m = 1.0. 





q — U 




N uni6riCcLl 


Semianalytical 


p 


A 






Re(tJo) 


-Im(t<Jo) 


1 
1 


iU 


z.you X iu 


y. ( oo X iu 


2.911x10-1 


9.780x10-2 




1 0-4 
iU 


z.yzo A iu 


y . ( o-i X iu 


2.910x10-1 


9.797x10-2 




iU 


9 Qi /I V 1 


y . ( zo X lu 


2.896x10-1 


9.771x10-2 




iU 


9 77n V 1 

Z. / / U X lU 


y.4oo X lu 


2.753x10-1 


9.490x10-2 






fi 1 c;q V 1 n~-^ 

o. -Luy A lU 


Q 1 9Q V 1 0^-^ 


8.144x10-2 


3.137x10-2 


9 


1 

lU 


A ftzin V 1 n~i 
^.o'y^u X lu 


y .ooo X iu 


4.832x10-1 


9.680x10-2 




1 n-4 

iU 


•4^.(5oo X lU 


o.y^o X iu 


4.830x10-1 


9.677x10-2 




iU 


4.(510 X iU 


o QQQ \/ 1 n— 2 
o.yyo X iu 


4.809x10-1 


9.643x10-2 




10^^ 


4.598x10 


8.880x10 


4.592x10-1 


9.290x10-2 




10-1 


1.466x10-1 


3.068x10-2 


1.466x10-1 


3.070x10-2 


3 


10"^ 


6.769x10-1 


8.662x10-2 


6.752x10-1 


9.651x10-2 




10-* 


6.754x10-1 


8.654x10-2 


6.749x10-1 


9.647x10-2 




10"^ 


6.732x10-1 


8.660x10-2 


6.720x10-1 


9.611x10-2 




10"^ 


6.437x10-2 


9.200x10-2 


6.428x10-2 


9.235x10-2 




10-1 


2.091x10-2 


3.054x10-2 


2.091x10-2 


3.056x10-2 



Table II: Fundamental frequencies for the electromagnetic field in SdS, obtained using the numerical and semianalytical methods. 
In this table, m = 1.0. 





q= 




Numerical 


Semianalytical 


I 


A 


Re(tJo) 


-Imiujo) 


Re(LJo) 


-lm(tjo) 


1 


lo-'^ 


2.481x10-1 


9.226x10-2 


2.459x10-1 


9.310x10-2 




10-* 


2.481x10-1 


9.223x10-2 


2.457x10-1 


9.307x10-2 




10-^ 


2.475x10-1 


9.176x10-2 


2.448x10-1 


9.270x10-2 




10-2 


2.374x10-1 


8.839x10-2 


2.352x10-1 


8.896x10-2 




10-1 


8.035x10-2 


3.027x10-2 


8.023x10-2 


3.033x10-2 


2 


10"^ 


4.577x10-1 


8.985x10-2 


4.571x10-1 


9.506x10-2 




10-* 


4.575x10-1 


8.991x10-2 


4.569x10-1 


9.502x10-2 




10-^ 


4.559x10-1 


9.439x10-2 


4.551x10-1 


9.464x10-2 




10-2 


4.371x10-1 


8.941x10-2 


4.364x10-1 


9.074x10-2 




10-1 


1.458x10-1 


3.037x10-2 


1.458x10-1 


3.038x10-2 


3 


10-^ 


6.578x10-1 


8.365x10-2 


6.567x10-1 


9.563x10-2 




10-* 


6.576x10-1 


8.349x10-2 


6.564x10-1 


9.559x10-2 




10-^ 


6.547x10-1 


8.399x10-2 


6.538x10-1 


9.520x10-2 




10-2 


6.276x10-2 


8.852x10-2 


6.267x10-1 


9.125x10-2 




10-1 


2.085x10-2 


3.039x10-^ 


2.085x10-1 


3.040x10-2 



The numerical simulations developed in the present 
work reveal an interesting transition between oscillatory 
modes and exponentially decaying modes. As shown in 
Fig. ^1 as the cosmological constant increases, the ab- 
solute value of — Im(a;o) decreases. 

Above a certain critical value of A we do not observe 
the exponential tail, since the coefficient k^xp is larger 
than — Im(wo) and thus the decaying quasinormal mode 
dominates. But for A smaller than this critical value, 
— Im(ci;o) turns out to be larger than k^xp, and the ex- 



ponential tail dominates. Certainly, for a small enough 
cosmological constant, the exponential tail dominates in 
the various cases considered here. 

Another aspect worth mentioning in the intermediate 
region is the dependence of the parameters fc|^p, 
fcg^p, and fcf°p on ^ and k^. The results suggest that 
the kexp are at least second differentiable functions of Kc. 
Therefore, close to Kc — 0, we approximate 

^:Sp(«c)«^(«c + c^=«?) , (59) 
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Table III: Fundamental frequencies for the axial and polar gravitational fields in SdS, obtained using the numerical and 
semianalytical methods. In this table, m = 1.0. 





q — U 




N uni6riCcLl 


Semianalytical 


p 


A 




-iirn^LJo ) 


Re(tJo) 


-Im(t<Jo) 


o 
z 


iU 


0. ( 00 X lU 


0.000 X lU 


3.731x10"! 


8.921x10"^ 




1 0-4 
iU 


0. / ( X iU 


o.ooU X iU 


3.730x10-! 


8.918x10"^ 




iU 


^ 791 V 1 


o.oOU X lu 


3.715x10-! 


8.888x10"^ 




iU 


O.ODD X lU 


o.Ooo X iU 


3.560x10-! 


8.572x10"^ 






1 17Qxin~^ 


Q 090 V 1 0^-^ 


1.179x10-! 


3.023x10-^ 


Q 
O 


1 

lU 


^ QQQ V 1 n~l 


Si ^i77\^ 1 0~2 
0.0 / / X lU 


5.992x10"! 


9.272x10"^ 




1 0-4 
iU 


o.yyo X lu 


S fi7fi V 1 0~2 
0.0 / X iU 


5.990x10"! 


9.269x10"^ 




1 

iU 


c Q70 V 1 
o.y / z X iu 


o.y /IX iu 


0.9DDX iU 


y.2o4X iu 




10"^ 


5.725x10"! 


8.695 xl0"2 


5.718x10"! 


8.874x10-^ 




10-1 


1.900x10-1 


3.030x10"^ 


1.900x10"^ 


3.032x10-^ 


4 


10-^ 


8.106x10-1 


8.810x10"^ 


8.091x10-! 


9.417x10-^ 




10-* 


8.102x10"! 


8.781 xl0"2 


8.087x10-! 


9.413x10"^ 




10"^ 


8.070x10"! 


8.799 xl0"2 


8.055x10"! 


9.376x10"^ 




10"^ 


7.733x10"^ 


8.714x10"^ 


7.720x10"! 


9.000x10"^ 




10-1 


2.564x10"^ 


3.034x10"^ 


2.563x10"! 


3.036x10"^ 



Scalar{/=1) Electromagnetic (/= 1) Axial and Polar Grav. (/ = 2) Scalar(/=1) Electromagnetic (/= 1 ) Axial e Polar Grav. {/ = 2) 




Figure 9: Exponential tail for the scalar and electromagnetic 
fields with 1=1, and for the axial and polar gravitational 
field with 1 = 2. In the graphs, m = 1.0. 



Figure 10: Approaching of — Im(ajo) and kexp, for the scalar, 
electromagnetic and gravitational fields, in the SdS geometry. 
Above a certain critical value of A (roughly 1.7 x 10"^, 4.0 x 
10"^ and 1.2 x 10"^ respectively, for the parameters used in 
the graphs), a tail it is not observed. For all curves, the mass 
is set to m = 1.0. 



(60) 

Previous results are illustrated in Fig. [TTl 



B. Reissner— Nordstrom— de Sitter black hole 

We assess here the behavior of the fields in RNdS ex- 
terior geometries that are not near extreme, nor close to 
the asymptotically flat limit. Direct numerical simula- 
tions and semianalytical (WKB) methods were largely 



employed to characterize the fields in this region. 

For scalar perturbations with £ = 0, the effective po- 
tential is not positive definite. As already shown in 
[Hill, solutions of Eq. ^ with £^0 lead to a con- 
stant tail. It was observed that in the SdS geometry there 
is a dependence of the asymptotic value on the ipilO^x) 
initial condition, in the context of a Cauchy type initial 
value problem. We have checked that the introduction of 
electric charge does not alter this picture. 

If £ > 0, we introduce the fields, and for £ > 1 
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Figure 11: Dependence of k^^p with Kc and I, in the SdS 
geometry. The symbols indicate the numerical values, and 
the solid lines are the appropriate fittings. For the left graph: 
Kip = 1.077 X lO""* + 0.984kc + 3.545k?, Klj, = -1.863 x 
10"* + 2.010kc + 2.608k? and kH^ = -1.959x 10"*+3.028kc + 
2.978«;^ For the center graph: fef^p = 2.082 x 10"'' + 1.988kc + 
6.141k?, klip = 2.712 X 10"'' + 2.974^^ + 8.284k? and kli^ = 
3.737 X 10"* + 3.984kc +4.517K?. For the right graph: fcjfp = 
2.616 X 10~* + 2.974kc + 9.895k? and fe"fp = 4.484 x 10"'' + 
3.896kc + 18.92k?. In the graphs, m = 1.0. 




200 400 600 800 1000 1200 1400 1600 1800 2000 



Figure 12: Quasinormal modes for RNdS and Z2 fields. 
The parameters for the geometry are q = 0.5 and m = 1.0. 
We have used 1 = 2. The results are similar for the other 
fields considered. 



the fields can also be analyzed. The first point stud- 
ied is the quasinormal phase. If the cosmological con- 
stant is high enough, the decay is dominated by quasi- 
normal modes, even when they no longer are accurately 
predicted by the expressions H53|) . This scenario, illus- 
trated in Fig. El is valid for all fields considered, with 
any charge smaller than its critical value. 

We have observed that the infiuence of the electric 
charge is mild, although not trivial. The range of varia- 
tion of the quasinormal modes with the charge is not very 



Figure 13: Dependence of Z^ field frequencies on q, for both 
the real and imaginary parts, in RNdS. The results are qual- 
itatively similar for the other fields considered. The parame- 
ters are m = 1.0 and A = 10~*. 



large. We illustrate this point in Fig. El H is interest- 
ing to compare the values obtained for the fundamental 
modes using the numerical and semianalytical methods. 
We find a very good agreement between these results. 
The difference is smaller for the first values of I. This 
is expected, since the numerical calculations work bet- 
ter in this region. In Tables ITVIVIII we illustrate these 
observations for a few values of A and q. 

For all fields considered, with a small enough cosmolog- 
ical constant, there is a qualitative change in the behavior 
of all fields considered. The late-time decay is dominated 
by an exponential tail. Therefore in the RNdS geometry 
we have 



with t ^ 00 



with t ^ 00 



with t 



(61) 



(62) 



(63) 



for t sufficiently large. At the event and the cosmological 
horizons t is substituted by v and u, respectively. Fig ure 
1141 illustrates this point, which was noted in [ij, [IJ for 
scalar fields, and we have extended this consideration to 
coupled electromagnetic and gravitational fields. In the 
aforementioned figure we compare the exponential tails 
at the event horizon for exterior RNdS geometries. 

An important point is that, unlike the quasinormal 
mode frequencies, the exponential coefficients have shown 
no dependence on the black hole's electric charge, for all 
kinds of fields at hand. Close to Kc = 0, our results are 
compatible with the expressions 



,(kc) «£(kc + c'''=k?) 



tp(K,)«(^-Hl)(Kc + C^^K?) I 



1,2 



(64) 



(65) 
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Table IV: Fundamental frequencies for the Z^^ field in RNdS, with g = 0.5 and m = 1.0, using the numerical and semianalytical 
methods. 





q — u.o 




IN urn6ric£Ll 




c • 1-1 

o61Tll&IlcLlytlC&l 






l-< O I / . 1 1 




l-< at / . 1 1 


1 TTl 1 / J \ 

-iiiH^u^o ) 


1 

1 


1 


9 71 V 1 


Q f^i V 1 n"^ 


98 V 1 
o.Zo X lU 


y .oo X lu 




lU 


9 71 -v 1 


y.ou X lu 


9S-V 1 

o.Zo X lU 


y.oo X lu 




1 rj — 3 


9 7n V 1 


Q /17 V 1 n~-^ 


97 V 1 
O.Z ( X lU 


Q /1Q V 1 n^^ 
y.^y X lu 






9 (^n V 1 

Z.UU X lU 


Q 1 1 V 1 
y. 1 1 X lu 


1 Pi V 1 

O. lO X lU 


Q 1 /I V 1 n^^ 
y. 1^ X lu 




10^1 

-LU 


1.16x 10~^ 


4 0^^ V 1 

4:. UO X _LU 


1.40x 10~^ 


4 n^^ V 1 

4:. UO X _LU 


9 
A 


1 

lU 


/I Q/i V 1 


Q 79 V 1 
y. ( z X lu 


/I Q/i V 1 
'i.y'^ X lu 


Q 71 V 1 
y. ( i X lu 




ir)-4 


A QQ V 1 


Q 79 V 1 0"-^ 
y . ( z X ±u 


A QJ. V 1 
^.y^ X J.U 


Q 71 V 1 
y. ( 1 X ±u 






/I Q1 v 1 n^^ 


y.oo X lu 


4 Q9 V 1 


Q R7 V 1 C\~^ 
y.D ( X lu 






4.73x 10 


9.31 X 10 


4.73X 10 


9.30x10 




10-1 


2.09x10"^ 


4.10x10-2 


2.09x10"^ 


4.16x10-2 


3 


10"^ 


7.05x10"^ 


8.95x10-2 


7.05x10"^ 


9.77x10-2 




10-* 


7.05x10^^ 


8.94x10-2 


7.05x10-^ 


9.76x10-2 




10-^ 


7.02x10"^ 


8.87x10-2 


7.02x10"^ 


9.73x10-2 




10-2 


6.77x10"^ 


8.50x10-2 


6.76x10-1 


9.35x10-2 




10-1 


2.98x10"^ 


4.10x10-2 


2.98x10-1 


4.10x10-2 



Table V: Fundamental frequencies for the field in RNdS, with q — 0.5 and m = 1.0, using the numerical and semianalytical 
methods. 





q = 0.5 




Numerical 




Semianalytical 


e 


A 




-Im(ci;o) 


Re(ajo) 


-lm(LJo) 


2 


10-^ 


3.82x10-1 


8.96x10-2 


3.81x10-1 


8.98x10-2 




10-* 


3.82x10-1 


8.96x10-2 


3.81x10-1 


8.98x10-2 




10-^ 


3.80x10-1 


8.93x10-2 


3.80x10-1 


8.95x10-2 




10-2 


3.66x10-1 


8.67x10-2 


3.65x10-1 


8.67x10-2 




10-1 


1.60x10-1 


4.06x10-2 


1.60x10-1 


4.05x10-2 


3 


10-^ 


6.13x10-1 


8.59x10-2 


6.12x10-1 


9.33x10-2 




10-* 


6.13x10-1 


8.59x10-2 


6.12x10-1 


9.32x10-2 




10-^ 


6.10x10-1 


8.59x10-2 


6.10x10-1 


9.29x10-2 




10-2 


5.88x10-1 


8.50x10-2 


5.87x10-1 


8.97x10-2 




10-1 


2.58x10-1 


4.07x10-2 


2.59x10-1 


4.07x10-2 



for any q lower than its extreme value. The dynamics of 
the fields in de Sitter spacetimes is therefore very different 
from similar cases in anti-de Sitter geometries, in which 
for high charge there is an abrupt change in the fields' 
decay |SJ. We have also explicitly assessed the behavior 
of the fields as q 0, comparing their quasinor- 
mal frequencies and exponential tails to those observed 
in the SdS fields. As anticipated, we found that the 
Z2 field behaves like the SdS Z~^ field and that the Z2 
field behaves like the SdS Z~ field, if the charge is small 
enough. We have observed that the RNdS fields Z^ tend 
smoothly to the SdS fields Z^. 



VI. APPROACHING THE ASYMPTOTICALLY 
FLAT GEOMETRY 

Scalar fields in the SdS geometry near the asymptoti- 
cally flat limit were studied in [ij, ■ f fii^ '^^^^ there 
is a clear separation between the event and the cosmo- 
logical horizons, such that 

6 = !S^I± > 50 . (66) 

A new qualitative change occurs in this regime, namely, 
a decaying phase with a power law behavior. Such a 
phase occurs between the quasinormal mode decay and 
the exponential decay phases. The field cannot be simply 
described by a superposition of the various modes, which 
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Table VI: Fundamental frequencies for the Zf field 


in RNdS, with q — 0.5 and m 


— 1.0, using the numerical and 


semianalytical 


methods. 










q = 0.5 




Numerical 


Semianalytical 


(. A 


Re(ajo) 


-Im(a;o) 


Re(ajo) 


-lm(LJo) 


1 10^^ 


2.71x10 


9.51x10-2 


2.69x10-1 


9.55x10-2 


10 


2.71x10-1 


9.50x10-2 


2.68x10-1 


9.55x10-2 


10"^ 


2.70x10 


9.47x10-2 


2.67x10-1 


9.51x10-1 


10~2 


2.60x10 


9.11x10-2 


2.57x10-1 


9.15x10-2 


1 n — 1 
10 


1.16x 10-1 


4.08x10-2 


1.16x10-1 


4.08x10-2 


2 10~^ 


4.94x10 


9.71x10-2 


4.93x10-1 


9.72x10-2 


10 


4.94x10 


9.71x10-2 


4.93x10-1 


9.72x10-2 


10~^ 


4.92x10 


9.67x10-2 


4.91x10-1 


9.68x10-2 


10~2 


4.73x10 


9.30x10-2 


4.73x10-1 


9.31x10-2 


10-1 


2.09x10-1 


4.16x10-2 


2.09x10-1 


4.10x10-2 


3 10"^ 


7.05x10-1 


8.95x10-2 


7.05x10-1 


9.77x10-2 


10-* 


7.05x10-1 


8.94x10-2 


7.05x10-1 


9.76x10-2 


10-^ 


7.02x10-1 


8.86x10-2 


7.02x10-1 


9.73x10-2 


10-2 


6.77x10-1 


8.50x10-2 


6.76x10-1 


9.35x10-2 


10-1 


2.98x10-1 


4.10x10-2 


2.98x10-1 


4.10x10-2 



Table VII: Fundamental frequencies for the field in RNdS, with q = 0.5 and m — 1.0, using the numerical and semianalytical 
methods. 





q = 0.5 




Numerical 




Semianalytical 


e 


A 




-lm{Luo) 


Re(ajo) 


-Im(Lj()) 


2 


10-^ 


3.82x10-1 


8.96x10-2 


3.81x10-1 


8.99x10-2 




10-* 


3.82x10-1 


8.95x10-2 


3.81x10-1 


8.97x10-2 




10-^ 


3.80x10-1 


8.93x10-2 


3.79x10-1 


8.94x10-2 




10-2 


3.66x10-1 


8.64x10-2 


3.65x10-1 


8.66x10-2 




10-1 


1.60x10-1 


4.13x10-2 


1.60x10-1 


4.05x10-2 


3 


10-^ 


6.13x10-1 


8.56x10-2 


6.12x10-1 


9.33x10-2 




10-* 


6.13x10-1 


8.56x10-2 


6.12x10-1 


9.32x10-2 




10-^ 


6.10x10-1 


8.56x10-2 


6.10x10-1 


9.29x10-2 




10-2 


5.88x10-1 


8.49x10-2 


5.87x10-1 


8.97x10-2 




10-1 


2.58x10-1 


4.07x10-2 


2.58x10-1 


4.07x10-2 



would imply a domination of the power law phase. This 
is illustrated in Fig. El 

The situation for RNdS cases obeying Eq. is pre- 
sented in Fig. El As can be seen in this figure, we have a 
perfect power law tail developing for large v when A = 0, 
as expected. For the RNdS exterior geometry with low A 
values, this power law tail appears quite clearly between 
the quasinormal zone and the exponential tail. 

With such data, we can speak of three different regimes 
in the field dynamics when one approaches the asymptot- 
ically fiat limit: first, a quasinormal regime, with its char- 
acteristic damped oscillations, followed by an intermedi- 
ate regime for which the power law tail is visible, and a 
late-time region for which an exponential tail dominates. 
This qualitative picture is valid for all fields considered. 



VII. CONCLUSIONS 

We have identified three regimes, according to the 
value of A for the decay of the scalar, electromagnetic, 
and gravitational (or in RNdS) perturbations. Near 
the extreme limit (high A), we have analytic expressions 
for the effective potentials and the quasinormal frequen- 
cies. The decay is entirely dominated by the quasinormal 
modes (as in j^ll), that is, oscillatory decay characterized 
by a nonvanishing real part of the quasinormal frequency. 

In an intermediary parameter region (lower A), the 
wave functions have an important qualitative change, 
with the appearance of an exponential tail. This tail 
dominates the decay for large time. Near the asymptoti- 
cally fiat limit (A <C 1), we see an intermediary phase be- 
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Figure 14: Tails of the Z+ {(. = 1) and Z+ {(. = 2) fields 
in RNdS. The parameters for the geometry are A — 10"'* 
and m = 1.0. The results are similar for the other fields 
considered. 
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Figure 15: Approaching the asymptotically flat geometry in 
SdS. Straight lines in the log- log graphs indicate power law 
decay. In the graphs, m = 1.0. 

tween the quasinormal modes and the exponential tail — a 
region of power law decay. When A = 0, this region en- 
tirely dominates the late-time behavior. 

Finally, for scalar fields with £ = a constant decay 
mode appears, and its value depends on the ■0^(0, x) ini- 
tial condition. Fig. reveals the appearance of the con- 
stant value for large and its dependence on 'i/'f (0, x). 
The value of (pQ falls below 10~^ for '0^(0, a;) = 0. These 
results are compatible with the analytical predictions of 



|14| . The analytical characterization of these regions and 
the corresponding critical values of A are crucial to a bet- 
ter understanding of these qualitatively different regimes. 

A very important aspect of the field dynamics in the 
RNdS geometries is that the influence of the electric 
charge on the field behavior, in general, was shown to 
be quite restricted. In particular, we observed no depen- 
dence of the exponential tail coefficients with q. This 
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Figure 16: Scalar (£ = 1), and Z2 fields (£ = 2), approach- 
ing the asymptotically flat limit in RNdS. The parameters for 
the geometry are g = 0.5 and m = 1.0. The results are similar 
for the other fields considered. 



contrasts with previous results obtained in the anti-de 
Sitter case |1| , where the charge plays a fundamental role 
in the tails. A deep understanding of this fact depends on 
new analytical asymptotic results along the lines of the 
ones obtained in (14 | . These points are now under inves- 
tigation. Also, since the presence of a charge implies an 
internal structure similar to that of a rotating black hole, 
our results might be interpreted as a broader universality 
of the frequencies here obtained. This question deserves 
further study. 
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